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while  Reissner's  mixed  variational  principle  (1984,  1986),  applied  to  the 
incremental  boundary  value  problem  yields  consistent  mixture  constitutive 
relations.  In  order  to  assess  the  model  accuracy,  numerical  experiments 
were  conducted  for  static  and  dynamic  loads.  The  prediction  of  the  model 
in  the  time  domain  was  obtained  by  an  explicit  finite  element  code.  DYNA2D 
is  used  to  furnish  numerically  exact  data  for  the  problems  by  discretizing 
the  deLails  of  the  microstructure.  On  the  other  hand,  the  model  capability 
of  predicting  effective  tangent  moduli  was  tested  by  comparing  results  with 
NlKE2d.  In  all  cases,  good  agreement  was  observed  between  the  predicted  and 
exact  data  for  plastic,  as  well  as,  elastic  responses. 


Abstract 

An  asymptotic  mixture  theory  of  fiber- reinforced  composites  with  periodic  microstructure  is 
presented  for  rate-independent  inelastic  responses,  such  as  elasioplastic  deformation.  Key  elements 
are  the  modeling  capability  of  simulating  critical  interaction  across  material  interfaces  and  the 
inclusion  of  the  kinetic  energy  of  microdisplacement.  The  construction  of  the  proposed  mixture 
model,  which  is  deterministic,  instead  of  phenomenological,  is  accomplished  by  resorting  to  a 
variational  approach.  The  principle  of  virtual  work  is  used  for  total  quantities  fo  derive  rrixtve 
equations  of  motion  and  boundary  conditions,  while  Reissner's  mixed  variational  principle  (1984, 
1986),  applied  to  the  incremental  boundary  value  problem  yields  consistent  mixture  constitutive 
relations.  In  order  to  assess  the  model  accuracy,  numerical  experiments  were  conducted  for  static 
and  dynamic  loads.  The  prediction  of  the  model  in  the  time  domain  was  obtained  by  an  explicit 
finite  element  code.  DYNA2D  is  used  to  furnish  numerically  exact  data  for  the  problems  by 
discretizing  the  details  of  the  microstructure.  On  the  other  hand,  the  model  capability  of  predicting 
effective  tangent  moduli  was  tested  by  comparing  results  with  NIKE2D.  In  all  cases,  good 
agreement  was  observed  between  the  predicted  and  exact  data  for  plastic,  as  well  as,  elastic 


responses. 
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Introduction 

The  increasing  use  of  fiber-reinforced  composites  in  the  engineering  technology  requires  the 
development  of  efficient  models  in  order  to  properly  simulate  their  stadc  and  dynamic  behavior , 
such  as  effective  moduli  and  wave  dispersion  effects.  Furthermore,  other  applications  may  also 
involve  load  amplitudes  that  induce  global  material  nonlinearities.  Consequently,  a  nonlinear 
material  model  is  needed. 

The  effective  stiffness  theory,  which  is  a  high-order  continuum  model,  was  developed  by 
Achenbach  and  Herrmann  (1968)  to  simulate  elastic  wave  dispersion  in  fiber-reinforced 
composites.  This  theory  has  been  later  applied  by  Bartholomew  and  Torvick  (1972),  Hlavacek 
(1975),  Achenbach  (1976),  and  Aboudi  (1981).  By  modifying  the  original  methodology,  Aboudi 
(1982,  1985)  extended  the  linear  moael  to  account  for  visco-plastic  material  responses.  In 
addition  to  the  effective  stiffness  theory,  mixture  models  have  been  developed  by  many 
investigators  (Martin,  Bedford  and  Stem,  1971;  Choi  and  Bedford,  1973;  Hegemier,  Gurtman 
and  Nayfeh,  1973;  Hegemier  and  Gunman,  1974;  Nayfeh,  1977;  Murakami,  Maewal  and 
Hegemier,  1979;  Nayfeh,  Crane  and  Hoppe,  1984;  and  Murakami  and  Hegemier,  1986).  To  date, 
high-order  mixture  models  have  not  been  extended  to  include  nonlinear  material  responses  for 
arbitrary  wave  motion. 

In  a  highly  heterogeneous  medium,  the  large  number  of  material  interfaces  renders  the  direct 
solution  of  the  problem  extremely  complicated  and  time  consuming.  In  order  to  alleviate  such 
difficulties  by  deriving  a  set  of  partial  differential  equations  with  constant  coefficients,  two  levels 
of  homogenization  methods  have  been  introduced.  They  are:  (1)  0(1)  homogenization  methods 
(Babuska,  1976;  Bensoussan,  Lions,  and  Papanicolaou,  1978;  Sanchez-Palencia,  1980;  Tartar, 

1977),  and  (2)  an  asymptotic  mixture  model  based  upon  an  0(e)  homogenization  method 
(Murakami,  Maewal,  and  Hegemier,  1981;  Murakami,  1985;  Murakami  and  Hegemier,  1986; 
Toledano  and  Murakami,  1987;  Murakami  and  Toledano,  1988).  0(1)  methods,  which  yield 
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effective  moduli,  are  non-dispersive.  Therefore,  simulation  of  wave  dispersion  requires  an  0(e) 

(at  least)  homogenization  model,  where  e  denotes  a  typical  ratio  of  micro-to-macrodimension. 

A  high-order  mixture  homogenization  model  of  fiber-reinforced  composites  with  periodic 
microstructure  is  presented  for  rate- independent  inelastic  responses,  such  as  elastoplastic 
deformation.  The  theory  is  based  upon  an  asymptotic  expansions  with  multiple  scales.  Key 
elements  are  the  modeling  capability  of  simulating  critical  material  interaction  and  the  inclusion  of 
the  kinetic  energy  of  microdisplacements. 

A  variational  approach  is  adopted  in  order  to  construct  a  mixture  model,  which  is 
deterministic,  instead  of  phenomenological.  The  principle  of  virtual  work  is  used  for  toic: 
quantities  to  derive  mixrure  equations  of  motion,  while  Reissner's  mixed  variational  principles  ( 
1984,  1986)  applied  to  rate  boundary  value  problems  are  used  to  yield  mixture  constitutive 
relations.  This  is  achieved  oy  supplying  the  variational  equations  with  appropriate  trial 
displacement  and  transverse  stress  fields  which  must  satisiy  certain  periMicity  and  normalization 
conditions.  The  accuracy  of  the  resulting  model  clearly  depends  on  the  accuracy  of  the  trial 
displacement  and  transverse  stress  fields,  which  may  be  found  accurately  as  solutions  of  a  system 
of  rate-boundary  value  problems  defined  over  a  unit  cell.  This  procedure  has  been  successfully 
applied  by  Murakami  and  Toledano  (1988)  to  predict  transient  responses  of  bi-laminated 
composites  and  is  extended  here  to  fiber-reinforced  composites. 

In  order  to  assess  the  model  accuracy,  numerical  experiments  were  conducted  for  static  and 
dynamic  loads.  The  prediction  of  the  model  in  the  time  domain  was  obtained  by  an  explicit  finite 
element  code.  DYNA2D  (Hallquist,  1982)  is  used  to  furnish  numerically  exact  data  for  the 
problems  by  discretizing  the  details  of  the  microstructure.  On  the  other  hand,  the  model  capability 
of  predicting  effective  tangent  moduli  was  tested  by  comparing  results  with  NTKE2D  (Hallquist, 
1986).  In  ail  cases,  good  agreement  was  observed  between  the  predicted  and  exact  data  for 
plastic,  as  well  as,  elastic  responses. 
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Formulation 


Consider  a  body  occupying  a  bounded  domain  Q  of  R3  with  smooth  boundary  dCl,  which 
contains  fibers  of  circular  cross  section  periodically  distributed  within  the  matrix,  as  shown  in  Fig. 

1.  The  fibers  occupy  a  volume  and  the  matrix  a  volume  .  A  rectangular  Cartesian 
coordinate  system  x,  is  selected  such  that  X\  is  parallel  to  the  fibers.  A  periodic  microstructure  in  the 
macrocoordinate  plane  (  X2,  X3)  is  defined  by  a  cell,  as  shown  in  Fig.  2  for  a  hexagonal  array. 


The  following  notation:  (  )(“),  a=l,2  will  denote  quantities  associated  with  material  a. 
Unless  otherwise  specified,  the  usual  Cartesian  indiciai  notation  is  employed  where  Latin  indices 
range  from  1  to  3  and  Greek  indices  from  1  to  2.  Repeated  indices  imply  the  summation 

convention.  Also,  ( )j  and  (  )  t  are  used  to  denote  partial  differentiation  with  respect  to  Xj  and  time 
t,  respectively. 

With  the  help  of  the  foregoing  notation,  the  governing  equations  for  the  displacement  vector  u 


(a)  and  stress  tensor  associated  with  the  ath-constituent  are: 


(a)  Equations  of  motion 


cr(a)  . 
j  1  •  j 


+  f  =  p  u  , 


(a) 

i  ,  1 1 


_  (a)  _  (a) 

<7  .  .  =  O  . 

Ji  ‘  J 


in  G 


(a) 


where  is  a  constant  body  force  and  p  denotes  mass  density; 
(b)  Rate  constitutive  relations 


•  (a)  (a)  (eg) 

a .  =  u  ,  . 

1  1  IJ  it  1 


) 


in  am 

in  n<a) 


(1) 

(2) 
(3) 


where  C(a^eP^  is  the  tangent  modulus  tensor, 
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(c)  Interface  continuity  conditions 


d)  (2) 

u  .  =  u  . 

1  I 


a  a)  v  (1)  =  a (2)  v  (1) 
j  >  j  j  i  ) 


on  A 


(4) 


(d)  Appropriate  boundary  conditions  along 

(e)  Initial  conditions  on  u  (°9  and  u  t(a)  at  t=0. 

Conditions  (a)-(e)  define  a  well  posed  initial  boundary  value  problem  on  Cl.  However,  due  to 
the  large  number  of  heterogeneities  in  the  medium,  the  direct  solution  of  this  problem  represents  a 
formidable  task.  To  circumvent  this  difficulty,  the  method  of  homogenization  which  is  based  upon 
a  multiscale  representation  and  interior  domain  asymptotic  expansions  is  adopted  here.  This 
procedure  has  been  successfully  applied  by  Murakami  and  Toledano  (1989)  for  inelastic 
bi-laminated  composites  and  is  extended  in  the  present  work  for  fiber-reinforced  composites. 


Multivariable  Field  Representation 

The  material  properties  are  periodic  in  the  x-space,  in  which  the  periodicity  is  characterized  by 
the  cell.  This  representation  implies  that  displacement  and  stress  fields  will  vary  according  to  two 
basic  length  scales:  (i)  a  macro  length  characteristic  of  the  body  size,  and  (ii)  a  micro  length 
characteristic  of  the  cell  spatial  dimensions.  Further,  in  most  applications  of  interest,  the  microscale 
is  much  smaller  than  the  macroscale.  Therefore,  it  appears  convenient  to  use  an  asymptotic 
technique  based  on  the  two  scales  involved  in  the  problem.  To  this  end  new  independent 
microcoordinates  are  introduced,  such  that 

xA  =  Xi/e  ,  i=2,3  (5) 

where  e  is  a  typical  ratio  of  micro-to-macrodimension.  As  a  consequence,  all  field  variables  now 

become  functions  of  both  the  macrovariables  x  and  microvariables  x  * 

F(xi,t)=F(x.,x.*,t,e)  (6) 

Spatial  derivatives  of  such  a  function  F  then  take  the  form 
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F ,  i  (  x  j  ’ 1 )  =  F  *  i  (xrxk*.t,e)  +  iF*.  (x.  ,xk*,t,e) 


(7) 


where  ( )  is  used  to  denote  partial  differentiation  with  respect  to  x^*.  For  notational  convenience, 
F*  will  be  written  as  F  in  what  follows.  It  can  be  seen  that  macrovariables  x  account  for  slow 
variations  of  any  material  function,  while  microvariabies  x*  account  for  rapid  variations. 

Applying  (5)  and  (7)  to  eqns  (l)-(4)  yields  a  new  set  of  "synthesized"  governing  equations, 
given  by. 


(a)  Equations  of  motion 


„(«)  ^  1  _(o) 

M.j  g  j*.J* 

(a)  (a)  (a) 

+  =P  Ui,tt  m 

QW 

(8a) 

o(ct)  =  a(a> 

J  l  l  J 

in  a(a} 

(8b) 

(b)  Constitutive  relations 

o‘“’=c‘“V?  {  e  k  |  ( u  '^ )  +  i  e  *  t  (  u <a> ) } 

1  ‘  e 

in 

(9) 

s*  (u,0))=  i(uf>. 
r  1  '  >  2  '  X  .  1 

^  ■  'a)  \ 

1*  +  U  1  .  1c*  ^ 

in  a(a) 

(10) 

(c)  Interface  continuity  conditions 

0)  ,  (2) 

u  =  u . 

1  l 

a)  (i)  (2)  a) 

a  .  v  =  a  : .  v  on 

!  ‘  J  J  ‘  J 

AI 

(ID 

The  hexagonal  ceii  is  modeled  as  two  concentric  cylinder;),  as  shown  in  Fig.  2.  Micropolar 


coordinates  (  r,  9  )  are  introduced  such  that 


v 


tan  9  =  x  ?  /x  * 


(12) 


With  respect  to  these  microcoordinates,  eqn  (12),  a  cell  domain  consists  of  subdomains  A(!)  and 
A<=>  occupied  by  the  fiber  and  matrix  respectively,  as  shown  in  Fig.  2: 


A 0)  =  {  ( r,  9  )  |  0  <  r  <  .  0  <  9  <  2a  ) 

Ai7)=  (  ( T,Q)\yfn^  <r<  1 ,  0<  9  <  2a  )  (13) 
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where  n(a)  is  the  volume  fraction  and  satisfies 

n':,)+nC)-l  (14) 

Further,  in  the  cell  the  interface  is  defined  by  r=(n(1))1/’2  and  the  boundary  by  r=l.  In  terms  of  the 
polar  coordinates  (12),  the  x*-periodicity  condition  for  a  hexagonal  cell  with  the  concentric 

cylinders  approximation  is  expressed  as 

F(x.r.0,t)  =  F(x,r,0  +  n:.t)  at  r=l  (15) 

The  x*-periodicity  condition  implies  that  all  field  variables  will  assume  equal  values  on 
opposite  sides  of  the  cell  boundary.  Therefore,  it  is  only  necessary  to  consider  a  single  cell  in  order 
to  determine  the  dependence  of  these  field  variables  on  the  microcoordinates  x*.  This  assumption  is 
adopted  for  elastic  as  well  as  inelastic  deformations  under  the  premise  that  a  typical  length  of  both 
elastic  and  inelastic  regions  is  large  compared  to  the  cell  dimension.  This  generalization  is  based 
upon  the  observation  that,  in  an  inelastic  domain,  tangent  moduli  of  the  two  constituents  may  be 
well  approximated  by  a  locally  periodic  function  of  space  variables.  Since  the  tangent  moduli 
depend  also  on  the  current  stress  state,  exact  periodicity  conditions  for  tangent  moduli  might  not 
hold,  in  general.  However,  und'  r  the  premise  that  the  variation  of  the  tangent  moduli  with  respect  to 
the  macrovariables  is  slow,  inelastic  tangent  moduli  may  be  approximated  by  a  locally  periodic 
function  of  space.  Finally,  it  should  be  pointed  out  that  synthesized  displacement,  stress  and  strain 
tensors  are  all  continuous  with  respect  to  the  macrovariables  x,  while  retaining  discrete 
heterogeneous  Soueture  with  respect  to  the  microvariables  x*  in  each  cell. 

Variational  Formulation 

The  construction  of  the  present  mixture  model  is  accomplished  by  applying  the  principle  of 
virtual  work  to  the  boundary  value  problem  of  total  quantities,  and  Reissner's  mixed  variational 
principles  (1984,  1986)  to  the  incremental  boundary  value  problem.  This  last  operation  allows  one 
to  derive  consistent  model  constitutive  equations. 

The  principle  of  virtual  work  applied  to  the  fiber-reinforced  composite  yields 
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f  T  V  f  r  /■  c  (a)  T  1  ~  Ja)T  (a)  2  (a)T  ,  .  (a)  (a)  (a)  .  , 

JL2^J{(5e  +-^e*  )ff  -  (f  -p  u[t)}dx* 

ft  03 1 A  ■*>  e 

+  -  J  l5u'")T  -  5  u(1;T )  T*  dA  ]dXjdx,dx3 

% 

=  J  [  ^  j  5u(a)T  vT(a)  dx*  ]  dA  (16) 

3ft  T  1=1  A  w 

where  A(a)  denotes  the  x*-domain  occupied  by  constituent  a,  vT(a)  represents  the  traction  vector 

acting  on  the  boundary  B\'j  with  outward  normal  v  and  T*  is  the  interface  traction  vector.  For 

arbitrary  variations  of  ulCt),  with  (3)  and  (10)  as  definitions,  the  Euler-Lagrange  equations  of  (16) 
include  eqns  (8)  and  appropriate  boundary  conditions. 

The  above  variational  equation  (16)  for  total  quantities  yields  any  order  of  homogenized 
models  with  corresponding  equations  of  motion  by  substituting  therein  appropriate  trial 
displacements.  However,  for  high-order  homogenization  models  the  associated  rate-constitutive 
relations  are  obtained  either  by  the  application  of  Reissner's  mixed  variational  principle  (  Reissner, 
1984,  1986)  or  a  set  of  experiments.  In  order  to  use  Reissner's  mixed  variational  principle,  it  is 
necessary  to  rewrite  the  incremental  constitutive  relations  (9)  in  terms  of  in-plane  strains  and 
transverse  stresses: 

a?  =  E  ,  j  e?  +  (  E  ,  :]a;  (17a) 


(ISa.b.c) 


3 


et  =  [e,,  +e*: 


e  3  3  ^  e  3  3  -e23  +  “e-' 


* 

2  3 


2  e ,  ,  +  2  e  * 


3  l 


3  1 


2  e  ,  .  +  2  ■*  * 


1  2 


(18d) 


Also  [  lT  is  the  transpose  of  [  ],  subscripts  p  and  t  denote  in-plane  and  transverse  quantities 
respectively,  and  the  matrices  [  Eaj3  ]  are  functions  of  the  elements  of  Cyj^.  Reissner's  mixed 
variational  principle  applied  to  the  incremental  boundary  value  problem  defined  by  (8)-(l  1)  yields 
f  r  V  f  r  ?  •  to)T  ■  (a)  ~  •  (c) T  ■  (a)  -  ■  (a)T  (a)  •  (a) 

j  L  L  J  (6eP  a?  +  det  CTt  - 611  p  3 


n 


a=l  a) 
A 


+  8 a [a)T (  e;:a>+  [Ef!]T  e  (a)  -  [  E'^  lat(a)  )  ]  dx' 


■i  [  [  (  5  u (2)T  -  5ua)T)T*  +  5  f*T  (u(2)-  u(1))]dA]  d  x1  d  x,  d  x, 

j 


=  J  (  J  5u(a)Tvf(a)  dx*  }  dA  (19) 

on  T  0=1  A » 

where  a  (Q)  is  the  acceleration  vector.  The  strains  are  ail  computed  from  the  trial  displacement  field 

by  (3)  and  (10).  For  arbitrary  variations  of  u^a)  and  d(a\  the  Euler- Lagrange  equations  of  (19) 

yield  the  rate  constitutive  equations  for  e:(«)  in  (17b),  as  well  as  the  rate  equations  of  motion,  (8). 
the  rate  boundary  conditions,  (d),  and  the  rate  form  of  (11).  Eqns  (17a)  are  considered  to  be  the 


definitions  of  dp(c0. 


The  mixed  variational  equation  (19)  with  appropriate  trial  functions  u(a)  and  dt(a)  yields  the 
incremental  constitutive  relations  for  a  given  continuum  model. 


Asymptotic  Analysis 


9 


The  asymptotic  technique  starts  by  assuming  the  following  expansions  for  the  displacements, 
stresses  and  tangent  moduli  (Lene  and  Leguillon,  1982): 

y  icr\ 

(20a) 

» 


u (a)  (  x  ,  x*,  t ,  e )  =  uJqJ  (  x  ,  x*,  t )  +  e  u  :“J  (  x  ,  x\  t )  +  e  2  u  (  x ,  x*.  t )  + 


0(a)(  x  ,  x\c ,  £)  =  icr ^  (x  ,x*,t)  (  x  ,  x*.  t )  +  ecr  (  x  ,  x*.  t)  +  ...  (20b) 


r  (a)  (ep)  r  (a)  (qp)  r  (a)  (ep)  2  r  (a)  (ep) 

^  (0>  EC(1)  E  (2)  '  - 


(20c) 


where  U(n)(2),  0^(2)  and  C(n)(2)  satisfy  the  x* -periodicity  condition.  In  the  sequel,  a  class  of 

inelastic  materials,  including  hardening  elastoplastic  materials,  which  admit  a  rate  potential  and  have 
positive  definite  tangent  modulus  tensor  is  considered.  Therefore,  in  the  expansion  (20c), 

C(0)(a)<ep)  is  also  assumed  to  be  symmetric  and  positive  definite,  which  obviously  holds  for  elastic 

responses  where  only  C(Q)(“>(eP)  =  C  (°)  is  required  in  the  expansion. 


Substituting  (20a, b,c)  into  (8)  -  (1 1)  and  grouping  terms  in  equal  powers  of  £,  a  sequence  of 
microboundary  value  problems  (MB VP's)  defined  over  a  unit  cell  is  obtained. 


MB  VP  for  0(£-2): 


_(o)  n 

or „  ■»  =  0 
l  >  (-0  ■  J* 


in  A 


(a) 


afa)  =C(a)(ep)  e*  (  (a). 

i  j(-l)  i  jkl(0)  ckl  ^  (0)  ; 


(i)  a) 
u  i  <p)  "  u  ■  (0) 


(i)  (i)  q  a) 

a  , ..  v  =  a ,,  v 
j  i  t-i)  j  ji  (-D  j 


in  A‘" 
on  A 


From  (21a,b)  the  operator  for  u^)(a)  may  be  expressed  as 

r  f  u  '  a) )  =  f  C  'a)  ^  e  *  (  u  ^ )  1  -  0 

^  u(0)  >  ~  1  ^  i  jk  i  p)  ek  i  (  U(0)  >  J .  j*  u 


in  A 


(a) 


(21a) 

(21b) 

(21c) 


(22) 


A  solution  of  (22)  is 


1  0 


(a)  ,  . 

U(0)  “  U(0)  ^  X  ’  :  ^ 


e*(u(0))  =  0 


(23a) 


_  W  n 


(23b) 


Equation  (23a)  states  that  u^q)  is  independent  of  x* 


MB  VP  for  CXe'1): 


_  («)  _  n 
aji(0).j*  0 


in  A 


(a) 


.(a) 


,  p(a)(ep>  u(a))+e*  (U(a))l 

°  i  j  (0)  ^  i  j  k  1  (0)  1  S  k  l  ^  U  CO)  ’  C  k  1  l  (1)  ;  J 


U(1)  =U(2) 

id)  id) 


(1) 


(1)  (2) 

a .  ,  ~  v  =  a : '  v . 

j  i  (0)  j  j  i  (P)  j 


(l) 


in  A 
on  A  j 


(a) 


Eqns  (24a, b)  imply  that 

(a) 
'(1) 


L(u(a),>  =  rc(a)(ep)  e*  ru(a)n  --C(a)(ep) 

L'\uns)-[  vi-kl(0)  ekl  (U(1))j..  Lij 


■  (  (a)  ^ 
j  It  1  (0)  .  j*  C  Ic  1  U  CD)  ' 


(24a) 
(24b) 
(24c  ,d) 

(25) 


Eqn  (25)  shows  that  uq)(“)  is  governed  by  the  same  operator  as  defined  in  (22),  except  for  the 
RHS,  which  vanishes  when  integrated  over  the  cell.  Therefore  the  integrability  condition  for  U(i)(a) 


is  satisfied.  The  form  of  the  forcing  term  in  (25)  suggests  the  following  expression  for  u(D(“> 


(a) 

id) 


(°>\  ~  P  4, 


(x,x*  t)  =  epq(u^)xP(x*) 


(26) 


where  XjN  is  x*-periodic.  The  substitution  of  (26)  into  (25)  yields  a  MB  VP  for  each  x^W  which  is 
continuous  over  the  cell.  These  problems  are  defined  up  to  a  constant  vector  with  respect  to  x*. 

Therefore,  it  is  convenient  to  choose  XjN  such  that  its  average  over  the  cell  vanishes 

Xfq(4)  =  £  J  X  f  q(  x*)  dx*  =  0  (27) 

d=l  ;o) 

A 


MB  VP  for  0(1): 


.«*) 


(a)  (a)  f  (a)  (a) 

1 t '  i  ‘  CT  j  i  (0) .  j 


ct  : ... .  .  =  p  u : ■  '  . .  -  f  -  ct  . 

j  I  (1) .  j*  K  I  (0)  . 


in  A 


(a) 


1  1 


(28a) 


(ct)  =c(a)(ep)  f  e  (  u (a) )  +  e  *  ( u (a) )  1 

i  j  (1)  ij  k  U0) 1  cic  l  ^  U(I)  ;  ckl  '•U(2)'J 


+  C  (“) (ep)  r  e  (  u (a) )  +  e  *  ( u  ^ )  1 
^  ij  k  1  (1)  L  e  k  1  ^  U(0)  1  C  k  1  (-Ua)  ;  J 


in  A 


(a) 


(28b) 


,  0)  _  (2) 
i  (?)  i  (2) 


(1)  a)  (2  1) 

ct..  mv  =ct:  '  ,v: 

ji  (1)  J  J  i  0)  J 


on  A 


(28c.  d) 


The  0(1)  homogenization  method  consists  in  solving  the  MB  VP's  defined  by  eqns  (24)  and 
(26)  for  U(i)(«).  The  0(1)  equations  of  motion  are  obtained  by  posing  the  integrability  condition 
on  (28a).  According  to  the  Fredholm  alternative  theorem,  the  problem  defined  in  (28)  has  a 
unique  solution  up  to  a  constant  vector  with  respect  to  x*,  if  the  operator  for  U(2)(a)  in  (28) 


satisfies  the  integrability  condition,  i.e  the  range  of  L  (  U(2)(a> )  is  orthogonal  to  its  kernel 


=  11(0)  (  x  - £ )  (f°r  example  see  Mardsen  and  Hughes,  1983;  Toledano,  1989). 

The  0(1)  equations  of  motion  can  also  be  obtained  by  substituting  the  trial  displacement 
function 


(a)  , 

U  (  X 


,t,e)  =  ui(0)(x  t)  +  £s  (x  t)xfq(xk.*) 


(29) 


in  eqn  (16)  and  keeping  only  (1)  terms,  which  implies  that  the  kinetic  energy  due  to  the  0(e) 
displacement  is  not  included.  In  order  to  improve  the  model  capability  in  predicting  dynamic 
responses  of  composites,  it  is  natural  to  include  the  kinetic  energy  due  to  U(i)(«).  In  addition,  the 


0(e)  homogenization  method  requires  the  knowledge  of  the  0(e)  transverse  stresses  by  solving 


the  MB  VP  defined  in  (28)  for  i^a),  and  then  computing  from  (28b).  However,  due  to  the 


complexity  of  the  MB  VP's  an  attempt  is  made  to  construct  a  reasonable  approximation  of 

without  solving  for  U2(°).  The  construction  of  the  present  0(e)  model  is  achieved  by  supplying 
the  variational  equations  (16)  and  (19)  with  appropriate  trial  displacement  and  transverse  stress 
fields  which  must  satisfy  the  interface  continuity  condition  (11)  and  the  x*- periodicity  condition 
on  the  cell  boundary,  eqn  (15).  As  a  result,  only  approximate  solutions  of  the  MBVP’s  are 
necessary.  This  procedure  yields  mixture  equations  of  motion  and  consistent  rate  constitutive 
relations. 


Trial  Displacement  and  Transverse  Stress  Fields 

One  starts  by  defining  average  displacements  for  each  constituent  as: 

. .  (ct)  .  .  1  f  ,  (a)  2  (a)  .  , 

U  (>  +  £“a>+£  uc>  >d 


V  ^ 

V  v(a> 


(30) 


Eqn  (26)  shows  that  Uq)(«)  is  excited  by  Up(0),q-  Therefore,  the  mixture  formulation  becomes 
more  tractable  by  introducing  microdisplacement  variables  which  represent  u^o)tq  +  uq(0),p  .  such 


that 


S|”<(x't)siWu 


(a)v(1)dA 
p  q 


V  J  Pt 


v  0)  dA 
(1)  q 


A, 


s  =  s 

p  q  q  ? 


The  following  trial  displacement  field  may  be  used 

uiCa)(xj,xk*.t.e)  =  U;a)(xj,t)+eSpq(x,.t)xfq(xk*) 


(31a) 


(31b) 


(32) 


In  (32)  U(a)  is  the  0(1)  displacement  associated  with  each  constituent,  while  Spq  represents  the 


amplitude  of  the  0(e)  displacement  microstructure.  It  is  noted  that,  for  general  periodic 
composites,  equations  (32)  imply  that  the  resulting  model  includes  eleven  macrodisplacements 


and  corresponding  equations  of  motion.  Further,  XjPQ  depend,  in  general,  on  tangent  moduli  of 
the  constituents  and  must  therefore  be  evaluated  at  each  load  increment.  However,  good  accuracy 


is  achieved  by  assuming  to  depend  only  on  the  microcoordinates  x*.  Such  an  approximation 
for  the  displacement  microstructure,  allows  for  the  application  of  homogenization  methods,  even 


to  nonlinear  problems,  without  having  to  solve  MB  VP's  for  % ;P3. 

For  an  elastic  analysis,  the  exact  solution  of  (24)  for  has  been  obtained  hy  Murakami 
and  Hegemier  (’986).  However,  the  present  mixture  analysis  becomes  more  tractable  by  adopting 
an  approximate  form  of  the  exact  solution  for  u^)(ct).  The  exact  solution,  together  with  (32), 


indicates  that  good  accuracy  of  the  0(e)  displacement  microstructure  is  obtained  by  introducing 
five  amplitude  functions  S;j  (i=l- 3,  j-2,3),  and  expressing  the  microstructure  with  two  functions 


of  r  and  9.  Eqn  (32)  then  reduces  to 


u  [a)  ( x  ,  x* ,  t ,  e )  =  U  (  x ,  t )  +  e  ( S  .  ( x ,  t )  cos  9  +  S  .  3  (x ,  t )  sin  9  ]  g (a)  (r)  (33a) 


where 


g(aV)  -Oi 


a+1 


» 


( r  -  -— ) 
r 


S  2  3  "  S  3  2 


(33b) 

(33c) 


In  (33b)  5ap  is  the  Kronecker  delta.  The  functions  g^(r)  cos  9  and  g(°)(r)  sin  9  are  continuous 
over  the  cell  and  satisfy  the  x*-periodicity  condition  (15)  and  normalization  condition  (27).  The 
capability  of  the  trial  functions  (33a)  to  simulate  harmonic  wave  propagation  in  elastic 
fiber-reinforced  composites  has  been  demonstrated  by  Murakami  and  Hegemier  (1986).  Eqns 
(33)  are  used  here  to  model  nonlinear  responses. 


Substituting  (33)  into  (24b)  and  considering  the  0(e2)  difference  of  the  average  transverse 

stresses,  the  following  0(1)  trial  stress  field  is  obtained 

5 

=  *  ia)(*,t)  +  -T-  [T(9)]t  (x,t)  (34a) 

r 

where 


[T(  9  )  ]  = 


cos  29 

cos  29 

sin  29 

0 

0 

cos  29 

-  cos  29 

sin  29 

0 

0 

0 

sin  29 

0 

0 

0 

0 

0 

0 

cos  29 

sin  29 

0 

0 

0 

-  sin  29 

cos  29 

In  (34a),  x^  and  t  are  stress  variables  defined  analogously  to  (18a).  As  an  0(£)  trial  stress  field 
one  may  use  the  following  (Murakami  and  Hegemier,  1986): 

a  =fg(a)«(Q(9)]P(x,t)  (35a) 


where 


[Q  (  9  )  ]  = 


0 

3cos  9 

sin  9 

0 

cos  9 

3sin  9 

0 

sin  9 

cos  9 

2sin  9 

0 

0 

2cos  9 

0 

0 

In  (35a)  P  is  an  0(e)  stress  variable  defined  by 

P.(x,t).-Lf  „<■%<'>  <*4  f<3(a,mv",dA 
1  e V  J  1  '  }  V  J  J  >0)  j 


! 


Integrating  (8a)  over  A(a)  and  making  use  of  the  x* -periodicity  condition,  the  following  mixture 
equation  is  obtained 


(a)  (aa)  ,  ,  ,  a  +  [  (a)  .  (a)  (a)  (a)  (a) 

n  a . .  .  +  ( -1 )  P.rn  r .  =  n  p  U. 

J1,JV  1  1  r  l.CC 


(37) 


where  the  average  operation  is  defined  by 

F (aa)  (  x ,  t )  h  — 1 


(a) 


I 


F  ^  (  x ,  x*.  t )  dx* 


A 


(a) 


(38) 


As  a  result,  the  trial  transverse  stress  field  may  now  be  written  in  the  form 

g 

ot(a)(x,x*,t,E)  =  r  (a)(x,t)  +  -f-[T(0)]t  (x.t) 


+  |  g (a)  (r)  [  Q  (0)  ]  P  ( x  ,  t ) 


(39) 


Mixture  Equations 

Substituting  (33)  into  (16)  yields  the  mixture  equations  of  motion,  associated  boundary  and 
initial  conditions,  while  the  rate  form  of  (33)  and  (39)  into  (19)  yields  the  model  rate  constitutive 
relations.  These  equations  are: 

(a)  Equations  of  motion 


+  (-l)aM  P  +n(a)f<<‘>  =  r,(Va>US‘) 

A  .  .  +  —  (a^  -  a?a)  +  R,.)  =  IS.,  ,  i  =  1.2 

j  i  .  j  2  2  i  2  i  2 i '  i2,it 

£ 

*j..l  ♦  4r<*J? •  »-U 
£  ' 


1  ..3 

2 


(40) 

(41a) 

(41b) 

(41c) 


where 


eCVr.  ,M{j  )h  i-  j  a  f“}  g  (ct)  ( cos  9  ,  sin  9 )  d  x* 

ct-=  1  'aj 
A 

RM  =  -4r  f  -T  (<*?!  c°s29  +  (jfi  Sin29)dx*  ,  i=l,2 
nU)AJAfflr- 


R 


=  — 1 —  f  —  ( -  a  ^  cos  20  +•  o'  ®  sin  29 )  d  x*  ,  i  =1 , 3 

1  _  (2) ,  J  . :  3  1  2 1 


n  A  -a»  r 


R,,= 


1 


23  n®A>2r 


f  1  (  n  ® 

J  TT(a22 


4-  a  )  sin  29  d  x* 


Inf  h(Va)  .  h(1)  =  I  ,  h  (2)  - - - —  ( 2  +  n  ®  in  n  (1^ ) 

Z!  4  4n(2)  n(2) 


(42) 


(43a) 


(43b) 


(43c) 


(44) 


In  (43)  A  (=rc)  denotes  the  area  of  the  cell, 
(b)  Boundary  conditions 


where 


5  U (<x)  = 

I 

0  OR 

„(,V")v.=T(0|,> 

J  1  J  l 

(45) 

5  S  ,  =  0 

i  2 

CR 

^  .  v  .  =  T  2  ,  i  =  1,2 

) 1  i  i 

(46a) 

5S. . =0 

i  3 

OR 

.'ijiV1'.’  ■  i=1'3 

(46b) 

5  S  ,3  =  0 

OR 

(iJjj+jajjJVj-T’+Tj2 

(46c) 

T  (ap)  =  JL 

‘  ”  A 

J  vt'“)  di* 

(47a) 

A™ 

e  ( T  2 ,  T 3 )  s 

v  t  l  7 

4il’ 

'  T  ':<x)  g (a>  ( cos  9  ,  sin  9  )  d  x* 

(47b) 

a“  am 

(c)  Initial  conditions 


ur>,u<*  ,S,2.S,3.Si:.,,Si;„  a,  t  =  0 


(d)  Consdtudve  relations 


J  [Ef2]dx*  x(0)  +  5a2  J  -i  [E?>  j[T]dx*t  +||g(a)[E^][Q]dx*P 


2  (<*)  ,  ( *!• )  c?  s  ,  f  rcWiT  ; 

n  "(eMt  +  - — S)  +  J  [Ei'J  e 

*1  (a) 

A 


<*>  H  v* 


|  -L  [T]T[E®]dx*  xa)  +  J  -L(T]T[E®  ][T]dx*t 

A«f 

+  fj  -Vs(2)[T]T[E;ff)2][Q]dx*P  =*[V]S 

A3,r' 

+  J  ~  £TlX[Ef*  ]Tef)l  d  x* 

f£  J  g(a)[Q]T[E?l]dx*  t  (a)  +  -^-g(2)[Q]T[E?l]  [T]d  X*t. 


2 

+  £- 
16 


£  j  g(“i2(Q!T(E“]tQ]dx*P  =  i(l)  ®-  0 


“’a"* 


wnere 


(a) 


(a) 


rr  (a) 


(a) 


(a) 


(a) 


—  rrT‘“''  r’  rr ,  rrW  rrw  .  r ,  (a)  ttw  ,  n  w  •>' 

~  (ID  t  "  [  U  2  .  2  L  3  .  3  U  2  .  3  +  L  3  ,  2  U  3  .  1  +  U  1  ,  3  U  1  .  2  +  U  2.  1  1 

S  =■  [S,,  S,,  2S„  S„  S,,]T 


(a) 


(a)  ,T 


2  2  3  3 


23  1  3 


(52) 


[V]=  — 
n (1) 


-1/2  1/2 

-1  -1 

0  0 

0  0 

0  0 


0  0  0 

0  0  0 

-1/2  0  0 

0  1  0 

0  0-1 


[R,]  = 


2  2  0  0  0 
0  0  0  0  0 
0  0  0  0  0 


’  0 

0 

0 

0 

2n 

'o 

0 

0 

2 

o' 

[  R  2]  = 

3 

1 

0 

0 

0 

.  [R3J  = 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

m 

1 

3 

0 

0 

0 

(53) 


■c(“(  and  t  are  defined  analogously  to  (18a).  It  is  noted  that  for  hardening  materials,  [  H22^cx->]  on 

the  LHS  of  (49)-(51)  is  symmetric  and  positive  definite.  Specific  forms  of  [Eij(cO]  for  the_von 
Mises  yield  criterion  and  associated  flow  rule  with  isotropic  strain  hardening  are  given  in 

Appendix  A.  The  solution  of  (49)-(5 1)  yields  t  (-a\  i ,  and  P  in  terms  of  e^te),  $,  eu(a),  and 

0(2)- Cri).  For  elastic  constituents  it  is  possible  to  evaluate  the  above  relations  explicidy 
(Murakami  and  Hegemier,  1986).  However,  for  inelastic  responses  the  integrals  in  (49)-(51)  have 
to  be  evaluated  numerically  at  each  increment. 

Substituting  eqn  (39)  into  (40)  and  (41 )  yields 


1  9 


2  2  3  3 

M22  =  3hP:/4  ,  M33  =M2]  =  hP2/4  ,  M33=3hP3/4  , 

2  3  3  2 

M,  2  =  M3i  =hPj/2  ,  M22=M23  =  hP}/4  , 

2  3 

jij,*M12  =  0  ,  h  =  h  (1)  +  h  (2)  (54) 

R2i=ti2/a(1)  ,  R22  =  (f22/2  +  t  33)/n(1)  .  R  3  t  =  - 1  3  { In (1)  , 

R23  =  t23/2n(1)  ,  R33=(-t22/2  +  t33)/n(1)  (55) 


The  remaining  constitutive  equations  for  CTj  j(cta)  and  Mu  are  obtained  from  (17a),  (38),  (39)  and 


(42).  The  results  are 


_  (°0  •  (aa) 

n  7C  cj  j  j 


=  j  EfJ  e dx*  +  J [ E  j°^]  dx*  i 


(a) 


(a) 


5a2J  ~  [E(i2!_3(T]dxM  +  |  j  g(a)[Ejai[Q]dx*  P 


(56) 


(a) 


(a) 


£*  M  p  =  X  J  S  (a)fl  E  m  e  f  i  <***  +  X  f  g  (0>  n  (  E  1°;]  dx*  "  1 


(a) 


0=1  A(a) 


^  J 

Cl=l  ^  (a) 
1 


+  SaJ  4g<I)n(E®](Tldx*i+£.2j 
.  <2>  r  a=i  „ 


g  (a)2  n  [  E  3“2]  ( Q  ]  dx*  P 


(57) 


0=1  (a) 

A 


where 

M  «  [  M  j  ,  M  , ,  ]T  ,  n  =  [  cos  9  sin  9  ]T  (58) 

The  accuracy  of  the  model  in  predicting  elastic  harmonic  wave  propagation  has  been  amply 
demonstrated  by  Murakami  and  Hegemier  (1986).  Therefore,  the  accuracy  of  the  model  will  be 
assessed  by  considering  inelastic  deformations  under  monotonic  and  step  loading. 


Reduction  to  an  0(1)  Model 


20 


The  0(1)  homogenized  model  can  easily  be  obtained  from  the  mixture  model  by  introducing 
the  kinematic  constraints: 


l 


U  . 
1 


(59) 


and  taking  the  limit  as  £  — >  0.  With  (59),  eqns  (40)  and  (41)  reduce  to 

_  (m) 


.(2 a) 

2  i 


a3  i  * 


-  (m)  (m)  . , 

i-  t  .  =  p  L  . 

1  r  1  . 

t  t 

(60) 

i=l-3 

”  ST 

+ 

73 

II 

o 

i=l  3 

(61a. b) 

where 


a=! 


(62) 


The  rate  constitutive  relations  are  obtained  from  (49)-(51)  with  (59)  by  setting  e=0.  The  results 
obtained  from  (49),  (50)  and  (43)  can  be  used  to  eliminate  S,j  at  each  increment  using  (61).  This 
procedure  yields  the  incremental  effective  stress-strain  relations  given  in  Appendix  B. 


Numerical  Results 

The  accuracy  of  the  present  model  is  tested  by  comparing  results  under  monotonic  and  step 
loading  with  numerically  exact  data  obtained  from  NTKE2D  and  DYNA2D.  In  what  follows, 
material  1  is  referred  to  the  fiber  and  material  2  to  the  matrix. 

For  monotonic  loading,  the  present  0(1)  model  is  tested  and  resuits  compared  with 

MKE2D  (Hallquist,  1986).  A  single  transverse  strain  component  is  considered,  say  LA^.  For  the 
0(1)  model,  only  one  typical  cell  (Fig.  2)  need  be  analyzed,  while  the  details  of  the  microstructure 
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are  discretized  in  NIKE2D  according  to  the  mesh  shown  in  Fig.  3.  The  material  properties  are 
given  in  Table  1.  Fig.  4  shows  the  global  response  of  the  composite,  where  the  effective 

transverse  stress  is  plotted  as  a  function  of  the  effective  strain  Uo^.  Excellent  agreement  is 

observed  between  the  0(1)  modei  and  NTKE2D.  Fig.  4  shows  that  the  homogenized  tangent 
moduli  are  E22^  =  12  GPa  for  the  elastic  range  and  EW™)  =  8  GPa  for  the  plasdc  range.  Figs.  5 

and  6  show  the  stress-strain  curves  at  the  center  of  the  fiber  and  in  the  matrix  at  r=l,  6=0, 
respectively.  Again,  good  agreement  is  obtained. 

Next,  the  dynamic  analysis  consists  of  the  waverefiect  problem,  illustrated  in  Fig.  7  wherein 
the  same  material  properties  are  employed  (see  Table  1).  The  half  space  is  subjected  to  a  step 
pressure  loading  of  intensity  1.4  GPa.  The  finite  element  mesh  used  in  DYNA2D  is  identical  to 
that  shown  in  Fig.  3,  except  that  now  20  cells  are  considered.  The  results  of  the  present  mixture 
model  are  obtained  by  subdividing  each  cell  into  39  two-node,  linear  isoparametric  elements  for 

the  elastic  computation,  and  11  elements  for  the  elastoplastic  case.  The  value  of  e  is  100v2 

microns.  The  time  variations  of  the  normal  stress  tJi2^  are  computed  at  the  center  of  the  fiber, 

and  in  the  matrix  at  r=l,  0=0  ,  of  the  eleventh  cell  from  the  impact  surface  For  an  elastic 
analysis,  results  are  shown  in  Fig.  3a  for  the  fiber  and  Fig.  8b  for  the  matrix,  while  for  an 
elastoplastic  analysis  they  are  shown  in  Fig.  9a  and  Fig.  9b,  respectively,  it  can  be  seen  that  a 
reasonable  correlation  is  obtained  between  DYNA2D  and  the  present  mixture  model. 

In  order  to  further  assess  the  accuracy  of  the  present  mixture  model,  the  waveguide  problem 
illustrated  in  Fig.  10  is  considered.  Due  the  symmetry  it  is  only  necessary  to  analyze  a  single  cell. 
TV  material  properties  are  given  m  Table  2.  The  composite  is  subjected  to  a  step  pressure  loading 
of  4.826  x  105  Pa  <  70  psi  ).  The  results  of  the  mixture  model  are  obtained  by  using  1 10 

two-node,  linear  isoparametric  elements.  The  cell  dimension  £  is  1.95  mm.  The  rear  surface 
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velocities  on  a  specimen  6.33  mm  thick  are  compared  with  the  experimental  data  recorded  by 
Hegemier,  Gurtman  and  Nayfeh  (1973).  Fig.  1 1  shows  the  time  variation  of  the  averaged  velocity 
for  an  elasdc  analysis.  In  the  case  of  an  elastoplasric  analysis  the  results  of  the  present  mixture 
model  are  compared  with  the  numerical  data  presented  by  Hegemier  and  Gurtman  (1974).  The 

volume  fractions  are  nH)=0. 16  and  n(2)=0.84,  while  the  cell  dimension  is  e=4mm.  The  average 

velocity  histories  at  Xi=5mrn  in  the  fiber  and  matrix  are  shown  in  Fig.  12a  and  12b,  respectively. 
CRAM  represents  the  resuits  obtained  by  a  two-dimensional  finite  difference  code,  while  TTNC 
represents  the  predicted  values  by  Hegemier  and  Gurtman  (1974).  It  can  be  seen  that  good 
agreement  is  obtained  between  the  available  data  and  the  proposed  mixture  model. 

Conclusion 

A  high-order  mixture  homogenization  model  of  fiber-reinforced  composites  has  been 
developed  to  simulate  rate- independent  inelastic  responses,  such  as  elastoplastic  deformation. 
The  construction  of  the  mixture  model  was  carried  out  by  supplying  the  principle  of  virtual  work 
for  total  quantities  with  an  appropriate  trial  displacement  field,  and  Reissner's  mixed  variational 
principle  for  the  incremental  boundary  value  problem  with  an  appropriate  trial  transverse  stress 
field.  The  former  operation  yielded  equations  of  motion,  boundary  and  initial  conditions,  while 
the  latter  yielded  consistent  model  constitutive  relations.  It  was  found  that  eleven  equations  of 
motion  were  required  for  the  mixture  model.  Good  correlations  with  experimental  and  numerically 
exact  data  indicate  that  the  proposed  model  may  be  used  to  simulate  critical  material  interactions 
for  static  as  well  as  dynamic  responses. 
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APPENDIX  A  -  The  Definitions  of  [Ey]  in  eqn  (14) 

It  is  computationally  advantageous  to  rewrite  eqn  (9)  in  terms  of  the  eiastoplasdc  compliance 
matrix  D(a)(eP) 


g  (“)  _  D  (a)(e'P)  a  1°) 


(Al) 


For  a  von  Mises  yield  criterion  and  associated  flow  rule  with  linear  strain  hardening,  D(°0(eP)  may 
be  expressed  as 


n(<x)(ep)  n  (a)  9  (a)  (a)T 

U  =  L)  +■ - —  s  S 


(A2) 


H'  is  the  strain  hardening  parameter,  ji  is  the  shear  modulus  and  the  deviatoric  stresses  s  and 


effective  stress  a  are  respectively  defined  by 


s  =cr . .  -  -*•  6  . .  o,  . 
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Using  (A2),  one  finds 


1 


c  _ 
E1  1- 


D 


(o)(ep) 


.  [E®]  — 


l  i 


D 


(a)(ep) 
I  1 


(  D„  D  j  3  D 


D ,  .  D 


[  £  f  \  ]  =  d  'a)(ep)  -  (  D  1 1  P  1  i  )  (aXeP) 


14  15 


i,j  =  2-6 


i  i 


(A3) 

(A4) 


(ct)(ep) 


(A5) 


APPENDIX  B  -  Effective  stress-strain  relations  for  fiber-reinforced  composites. 
The  following  stress-strain  reladons  are  obtained: 
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where 


D  SI  =  £  "  (a)  E  f?  +  a (1)  [  E  [  E  [  E  f  J]T  +  a (2)  {  [  E  fj]  [  M ,]  [  E  fJ]T 
a=l 

+  [E(1=5][M2]  {a}  +  {a}T[M2]T[E^)]T+{a}T[M3]  {a}  )  -{b}T[D(n)]  {b}  (B2) 


r  ^  (m).  (1)  r  _  (la).  r  _  (ta).-l  (2)  f  r  _  (2a).  ...  .  ,  .T  .  .  ,  ,T  ) 

[  D  j ,]  =  n  [E12][E22]  +  nw  l  [E12][Mj]  +  {a}  [  M  ,]  j 


-  {b}T[DCn)]  [A  ] 


(B3) 


[D^J]  =  n (1)  [  E 2  2Y*  +  n (2)  [ M  j]  -  [  A  f  [ D  (n)]  [  A  ]  (B4) 

(b}  =  [  M  J  [  E  f?]T  +  [  M  ;]  {a}  -  [  E  E  f?]T 

*[V]T  {  [M2]T[E(12:)]T  +  [M3]{a}  )  (B5) 

[A^^ME^-tvftM/  (B6) 

[D(n)]=  {-j-[E^2)rl  +  -^-([M1]  +  [V]T[M3][V] 
n  '  nu 

-[M;][V]-[V]T[M2]T)  }  '  (B7) 

[M  2]  =  [  E^V1  (  [I]-[Bj][M2]T)  ,  [M2]  =  -[Ef2)rl[B1][M3J  , 

[M3]  =  ((B2]-[Bl]T[E^)rl[B1])'1  (B8) 
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[B,]=-4t  f  4[Ef’][T]dx*  ,  ta,]  =  -ip  |  -!j  [T]T[E®](T]dx*  , 
n  wtc  a  r"  n  ic  » r 

A  A 

4(TlT(Ef’2]Tdx*  (B9) 

n  rc  r“ 

In  eqn  (B8a),  [  I  ]  denotes  the  identity  matrix,  {a}  and  {b}  are  5x1  column  vectors,  [  and  [ 
Di2^m^]  are  1x5  row  vectors  and  D^N  is  a  scalar.  All  other  quantities  are  5x5  matrices.  In 
addition,  matrices  [  I>22(m)].  [  and  [  M3]  are  symmetric. 
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Table  1  Material  Properties  for  Wavereflect  Problem. 


Material 

Volume 

Fraction 

Density 

Young's 

Modulus 

Poisson's 

Ratio 

Yield 

Stress 

Hardening 

Parameter 

a 

n(“) 

p(ct)(kg/m3) 

EW(GPa) 

v(«) 

aY(GPa) 

H'(GPa) 

1 

0.25 

2200 

69 

0.17 

2 

0.75 

1266 

6.395 

0.25 

0.7 

0.35 

Table  2  Material  Properties  for  Waveguide  Problem. 


Material 

Volume 

Fraction 

Density 

Young's 

Modulus 

Poisson's 

Ratio 

Yield 

Stress 

a 

n(«) 

p(a)(kg/m3) 

E<°>(GPa) 

v(“) 

ay(GPa) 

1 

0.272 

1850 

29 

0.38 

2 

0.728 

1290 

8.25 

0.36 

0.133 

matrix 


I  A  fiber-reinforced  composite  and  the  coordinate  system. 
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Fig.  3a  Tine  variation  of  nornal  stress  r,  for  an 
elastic  analvsis 
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